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Fra le formale per le quadrature sono notevoli quella dei 
trapezii 



(«) 



ove 



ff{x) dfc^^-^ [f{a) + fib)] + B, 



12 ^"^«)' 



e quella di Simpson 

b 

a 

ove 



6 1' 


w- 


i-*i 


'\V 


B=- 


(h- 

' 4! 


-ay 

'5\ 


>»>'(«), 



rappresentando costantemente cpp u u^ yjiJorQ jiplerj^^^W ^^ ^ 
e 6 (*). Il resto nella (a) è nullo, se f{x) è intera di primo grado, 
e nella (|3) il resto è nullo, se /*(«) è intera di grado non su* 
periore al terzo. 



. (*) Pabl)licai queir espressione del resto dellf formoU di Sl)M^e^0N ^elle 
Apj^iccuioni geometriche del calcolo infinitesimale , pag. 206. 
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i 6. f^EANO 

Parallelamente a queste formule si hanno quelle dì Gausd. 
L'analoga alla (|3) è 



b 



+ f 
OTe 



180.4»/ ^ ^' 

e .il resto è nullo per le funzioni di grado non supcriore al terzo. 

Paragonando le formule (|3) e ((5'), che si possono conside- . 
rare come egualmente approssimate, risulta che è più semplice, 

in generale, il calcolo dei tre valori f{a), /*(— 5 — )» f(^)^ che 
esige la formula (|3), che il calcolo dei due 

(a + h 1 b^a\ /a + & 1 h^a\ 

n 2 v^ 2 r\ 2 "^yi 2 ; 

che esige la iormola (|3'). Questo spiega il maggior uso della 
formula ((S) di Simpson sulla corrispondente {^') di Gauss. 

Le formule di Gauss costituiscono una successione infinita , 
mentrechè le formule dei trapezii e di Simpson erano finora 
isolate. Io mi propongo di esporre qui una successione di infinite 
formule di quadrature, di cui le due prime sono appunto la (a) 
e la (^). 

Per semplicità supporremo i limiti dell'integrale eguali a — 1 
e + 1; poiché tasta fare il cambiamento 

a+b b—a . 
X = — ^r— H — ar , 



onde ridurci a questo caso. 

La questione che ci proponiamo è questa: Determinare gli 
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f 

w+ 1 coeflicienti -4o» A» • • • -4,., e gli n— 1 valori Xi x^ . . . 
Xj.^1 compresi fra — 1 e +1 in guisa che la formula 

(1) j jmdx = AJ{-l)+AJ(x,)+A,fix,) + . . . 

r ' +A^a(x^,) + Am 

sussista, qualunque sia la funzione f{x) intera di grado 2n — 1 . 
La soluzione è la seguente. Pongasi. 



(2) 



(rv-»"- 



Avendo la funzione (ir* — 1)" le radici — 1 ed 1 multiple 
d'ordine w, la sua derivata (n— !)"•", r„, avrà le radici a?o=— 1, 
a?„= 1 , semplici, ed w — 1 radici x^ x^ , , . x^^t distinte e com- 
prese fra — 1 e +1. Si calcolino i coefficienti J©, A."- colla 
formola 



+ 1 



— 1 

Allora sussisterà la formula (1). 

Infatti, si divida la f{x), funzione intera di grado 2n — 1, 
per F„, di grado w + 1 ; siano (f(x) il quoziente, ^(x) il resto, 
onde : 

(4) f{x) = ^{x)+Y„cp{x). 

Sarà (/* (x) di grado n, e cp (x) di grado n — 2. Attribuendo 
ad ic gli n+ 1 valori iPo, a?! . . .a?„, per cui si annulla F„, si avrà 

.Quindi la funzione ^{x), intera, di grado w, di cui si co- 
noscono i valori per n + 1 valori della variabile , si può espri- 
mere colla formola d'interpolazione di Lagrange: 

(5) •ut)=y (^-^o)- • • (^-^^*) (^^ - ^^i) • • • (^ - ^") / (^ ) 
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é G. t>EAKO 

Dalla (4) n ha: 

-ri -M +1 

(6) {f{x) da?=[^ (ir) dx+JT„(f{x)dx 
-1 ^1 -I 

Ora, dalle (5) e (3), li ncaira 

(7) j^{x)dx^2j,f{x^. 
—1 

Riguardo al secondo integrale, coirìntegrazione per parti si ha: 

(d \"~* / d \""^ 

=t f <"-»)(ic) . (ic»-- 1 )"iF fy ("-*)(ic) (x»— 1)" *; . 

Mettendo i limiti — 1 e + 1 , tutti i termini integrati nel 
secondo membro si annullano, perchè contengono il fattore a?* — 1; 
e siccome (f (x) e di grado n — 2, sarà (p(""'^(a?)==0, qx^^: 

■i-i 

(8) JY„f{x)dx=0 

Sostituendo nella (6) ai due integrali de} secondo meiìit^ro i 
joro valori dati dftUe (7) ed (8), si h^ la fonaola (1) obe bì 
voleva dimostrare. 

La formola (1), esatta se f{x) è intera di grado 2w — 1, 
è apprps9imata se f{x) è una funzione arbit^aida. Per aaloolare 
l'errore Mt tale che si abbia: 

+1 

(9) (fix)dx = lAJ{x;^+B, 



Digitized by 



Google 



GENERALIZZAZIONE DELLA FORMULA DI SIMPSON 7 

si formi la funzione F{x), intera, di grado 2n — 1, che Soddisfa 
alle 2n condizioni: 

Fia;,)^f{Xo), F{x,) = f{x,), F{x,) = f(x,), 

. . . F{x^,)=f{x^,) . J'(a?„) = /-(a:„), 

F' {X,) = f (X,) , F>{x,) =f'{x,),... F'{x^)=f'{x^,) , 

Si avrà, com'è noto: 

+1 

Integrando si avrà appunto \F{x)dx=^2^rf{^r)^ oi^de 

ri 

(11) B=\{x-Xo){x-x,y...{x-x^,y{x-x„)^-^dx. • 

—1 

Portando fuori del segno integrale il fattore ,^ ,, , cosa 

(2n) ! 

lecita, poiché il fattore rimanente ha un segno costante neirin- 

tervallo di integrazione, si ha : 

(12) ^=^r r^-'^o) (*-*•>' • • • (^-^-i)*(*-«-) *»• 

— 1 

Facendo n = l, si ha la formula dei trapezii (a). 
Per w = 2 si ha la formula di Simpson (|3). 
Per n = 3, fatti i calcoli, si ha: 



— 1 



3.5*. 7.6!' ^'' 
e il resto è nullo per le funzioni di grado inferiore al 6", 
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8 G. PEANO - GENERALIZZAZIONE DELLA FORMULA DI SIMPSON . 

Per w =: 4 si ha : 1 

I 






OTe 



2». 3*. 5*. ?• "" 2778300 



Nota. '— La stessa questione fa già trattata dal compianto D. Turazza, 
nel suo scritto : < Intorno alVuso dei compartimenti diseguali nella ricerca del 
valore numerico di un dato integrale > (Memorie L R. istituto Veneto, voi. V, 
^855, pag. 277-298). Credo però nuove le espressioni di Y^ e dei resti. 
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